Effect of a variation on the vibrational spectrum of Sr2 
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We consider the effect of a variation on the vibrational spectrum of Sr2 in the context of a planned 
experiment to test the stability of fi = me/rUp using optically trapped Sr2 molecules [Zelevinsky 
et ai, Phys. Rev. Lett. 100, 043201; Kotochigova et ai, Phys. Rev. A 79, 012504]. We find the 
prospective experiment to be 3 to 4 times less sensitive to fractional variation in a as it is to 
fractional variation in ^. Depending on the precision ultimately achieved by the experiment, this 
result may give justification for the neglect of a variation or, alternatively, may call for its explicit 
consideration in the interpretation of experimental results. 
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I. INTRODUCTION 

In the endeavor to understand nature on its most fun- 
damental level, physicists are striving for a description 
of all fundamental forces within a single unified theory. 
Some promising theories suggest that observable quanti- 
ties such as the electron-to-proton mass ratio /i = rrie/mp 
or the fine structure constant a = /he may not have 
fixed values A detected drift in /i or a could thus 
provide valuable insight into the fundamental workings 
of nature beyond our current understanding. To date, 
laboratory measurements have verified the stability of /i 
and a on the fractional level of 10~^^ [ll and 10~^^ Q per 
annum, respectively. More stringent laboratory tests are 
further motivated by recent evidence of a spatial gradient 
in the value of a based on an analysis of quasar absorp- 
tion spectra Q. It has been suggested that the Earth's 
motion relative this gradient may lead to measureable 
effects in the laboratory 

The Ye group at JILA (Boulder) aims to test the sta- 
bility of /X to high precision using Sr2 molecules confined 
in an optical lattice Q. The experimental protocol, as 
outlined by Zelevinsky et al. 0, Q , calls for optical Ra- 
man spectroscopy between select vibrational levels of the 
X ground electronic potential. A variation in /i al- 
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tcrs the vibrational spectrum, and experimental sensitiv- 
ity to this change may be optimized with a prudent choice 
of levels to incorporate into the spectroscopic scheme. It 
was shown in Refs. 0, that, with respect to variation 
in /i, the lowest and highest vibrational levels experience 
minimal displacement relative to the potential, whereas 
levels in the intermediate part of the spectrum experi- 
ence a much larger shift. Zelevinsky et al. have focused 
on transitions between the n = 27 intermediate level {n 
being the vibrational quantum number) and "anchor" 
levels at the bottom and top of the spectrum. Fig. [T] 
illustrates the basic objective of the experiment. A de- 
tected drift in the frequency ratio R (see Figure) is to 
be interpreted as a drift in the clectron-to-proton mass 
ratio. 

Here we investigate the effect of a variation on this 
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FIG. 1. (color online) Basic illustration of the experiment 
described by Zelevinsky et al. 0, 01 to test the stability of 
the electron-to-proton mass ratio using the vibrational spec- 
trum of **Sr2. The curve represents the X electronic po- 
tential, with horizontal dashed lines corresponding to select 
vibrational levels supported by this potential. Frequencies p 
and v' axe to be measured by an optical Ramsey scheme in- 
volving an intermediate excited electronic state (not shown). 
The dimensionless ratio R is sensitive to variations in ^ and 
is independent of any external reference (clock) frequency. 



promising experiment. The electronic potential depends 
on a through relativistic effects of electron motion. A 
variation of a alters the potential and, consequently, the 
vibrational spectrum supported by it. Thus, a measured 
drift in R may be due (or partially due) to a varia- 
tion, threatening misinterpretation of the experimental 
results. 

The experiment described by Zelevinsky et al. is sim- 
ilar in spirit to an experiment posed simultaneously by 
DeMille et al. Q to test the stability of ^ using diatomic 
molecules. A key difference is that, whereas Zelevinsky 
et al. focus on a single electronic potential, DeMille et 
al. suggest probing the splitting between vibrational lev- 
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els supported by different electronic potentials. Recently, 
we analyzed the influence of a variation on the experi- 
ment of DeMille et al. using a semi-classical (WKB) ap- 
proach, specifically focusing on the system Cs2 [lOl- We 
found the experiment to be order-of-magnitude as sen- 
sitive to fractional variation in a as it is to fractional 
variation in ^. Considering the anticipated precision of 
this experiment Q , together with the current laboratory 
limit on a variation [3|, we concluded that a variation 
may not be negligible for the proposed experiment. This 
finding largely motivated our present work. 



II. PRELIMINARY SET-UP 

Strictly speaking, only variations in dimensionless 
quantities have physical meaning. For a given vibrational 
level n, we will concern ourselves with the normalized en- 
ergy £n = En/ D, where En is the vibrational energy rela- 
tive to the bottom of the potential and D is the potential 
depth. Clearly, f„ is limited to the range < f „ < 1. 
Variations in fi and a induce a shift in £„, 
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where we employ the shorthand notation 
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The quantities d^£n and da£n quantify the sensitivity 
of the vibrational level n to fractional variations in the 
electron-to-proton mass ratio and the fine structure con- 
stant, respectively. 

We could, if desired, regard 5En — 5£n x _D as an 
"absolute" energy shift, an association which amounts 
to arbitrarily assuming the potential depth to be fixed 
with respect to any variation. Fixing any other energy 
reference — such as the atomic unit of energy, given by 
e^rrie/h^ = a^mgC^, or the SI unit of energy, which it- 
self references the hyperfine frequency of Cs as well as 
a platinum-iridium prototype mass held in Paris j7l| — 
would be equally justified and would generally yield a 
different "absolute" energy shift. Here we actively avoid 
the possibly slippery notion of absolute energy shift and 
quote results for and da£n, as these are unambigu- 
ously defined. 



III. THE MORSE POTENTIAL 

The Morse potential represents an idealized electronic 
potential for a diatomic molecule. It is given by 



V{r) = D 
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where r is the internuclear separation, with rp being the 
equilibrium distance, and is directly related to the 
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FIG. 2. The function }{x) = 2[x~l + VT^J . For the 
Morse potential, this function modulates the level sensitivi- 
ties to both ji and a variation across the vibrational spectrum, 
with the argument x taken as £n = En/D. The most deeply 
{£n — > 0) and loosely (f„ — >■ 1) bound levels are insensitive to 
variations, while intermediate levels have much larger sensi- 
tivities. 



width of the potential. The normalized vibrational en- 
ergies for the Morse potential are given precisely by the 
formula 



£n^ ^{t 



(1) 



where e = ha^j2/ DM and M is the reduced nuclear 



mass. 



From Eq. ([l} we see that a variation in may be at- 
tributed solely to a variation in the parameter e. Specif- 
ically, we may write 
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where f{x) = 2 [x — 1 + \/l — x\ . The function f{x) is 
displayed in Fig. [2] This function modulates the sensi- 
tivity of the various levels of the vibrational spectrum 
to variations in ji and a. Notably it approaches zero 
in the limits a; — >■ and x — > 1 and has a maximum 
at a; = 3/4. This translates to minimal sensitivities for 
the lowest and highest vibrational levels, with the largest 
sensitivities occurring for levels in the intermediate part 
of the spectrum. 

We may go a step further and, based on physical rea- 
soning, deduce a numerical value for the factor (9^ In e) 
appearing in Eq. This is accomplished most trans- 
parently by assuming atomic units, though we reiterate 
that e itself is dimensionless. When expressed in atomic 
units, the molecular potential (and its depth, width, 
etc.) is independent of the electron-to-proton mass ratio, 
whereas the reduced mass has a value which is inversely 
proportional to /i. From the definition of e, it follows 
that {d^Me) ^ 1/2. 
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In contrast to (9^1ne), there is not a simple analyti- 
cal result for the factor {da\ne). Nevertheless, we may 
provide a physically reasonable estimate of (i9alne) by 
realizing that, in atomic units, the electronic potential 
is independent of a in the nonrelativistic limit, hav- 
ing relativistic corrections which scale as (aZ)'^, with 
Z being the nuclear charge number. This suggests that 
(da Ine) ~ (aZ)'^. From this reasoning, we may suspect 
the vibrational spectrum of Sr2 to be nearly as sensitive 
to a variation as it is to ^ variation. 

The exercise of this section provides us with useful 
insight which is applicable to real physical systems. A 
true potential, of course, is not restricted to the form of 
a Morse potential. Nevertheless, the function f{x) dis- 
played in Fig. [2] is expected to give a qualitatively accu- 
rate depiction of the sensitivities and i9q£„ versus 
the normalized energy £„. In the vicinity of the equi- 
librium distance, the potential resembles that of a har- 
monic oscillator. The lower portion of the energy spec- 
trum is then well-described by a single term, proportional 
to (n -|- i) , in a Dunham-type expansion. Across this re- 
gion, d^£n vs £„ and daSn vs £„ arc essentially linear. 
Approaching the dissociation limit, anharmonic effects 
become important and the remaining terms in the Dun- 
ham expansion, proportional to {n+ ' + i) ' '^^'^■^ 
then drive the sensitivities back to zero. The fact that 
the sensitivities approach zero in the limits f „ — > and 
£n — > 1 is a consequence of our choice for the zero of 
energy (bottom of the potential) and our choice to nor- 
malize energy to the dissociation energy. 



IV. AB INITIO CALCULATIONS FOR LEVEL 
SENSITIVITIES 

We have calculated the A"^I]+ potential of Sr2 
using the relativistic computation chemistry program 
DIRACIO In order to reduce computational ef- 

fort we employed the infinite order two-component rel- 
ativistic Hamiltonian obtained after the Barysz-Sadlej- 
Snijders (BSS) transformation of the Dirac Hamiltonian 
in a finite basis set [isj . This approximation includes 
both scalar and spin-orbit relativistic effects to infinite 
order and is one of the most computer time efficient and 
accurate approximations to the four-component Dirac- 
Coulomb Hamiltonian. Electron correlation was taken 
into account using closed-shell single-reference coupled- 
cluster theory including single, double, and perturbative 
triple excitations [CCSD(T)]. The Faegri dual family ba- 
sis set [l^ was used, augmented by diffuse and high 
angular momentum exponents to obtain 21sl8pl2d6f2g 
Gaussian orbitals. Virtual orbitals with energies above 
45 a.u. were omitted, and the 56 outer core electrons were 
correlated. 

We subsequently fed our CCSD(T) potential curve into 
a Matlab routine to solve the Schrodinger equation for 
the nuclear part of the molecular wave function within 
the Born-Oppenheimer approach. A symmetric three- 



point finite difference method was applied to obtain the 
nuclear eigenfunctions together with their corresponding 
vibrational energies. For the discretization of the inter- 
nuclear distance a step size of 7 X 10-3 A was chosen. 

We may gauge the accuracy of our ab initio method 
through direct comparison with experimental results of 
Gerber et al. [l^. These authors have tabulated en- 
ergies for the n = through n = 35 portion of the 
X^S^ vibrational spectrum. Furthermore, they have 
determined the dissociation energy of this state to be 
D — 1060(30) cm"-'^. Our computed dissociation energy, 
D = 993 cm^^, is about 2a lower than the experimen- 
tal value. Comparing individual levels, we find that our 
computed vibrational energies differ from experimental 
values by no more than 2 cm"^ for levels spanning the 
lower half of the potential depth (n = through n ~ 15). 
Above this, our computed energies steadily diverge from 
experiment values, with our values being increasingly 
smaller in comparison. For example, for n = 21 our 
computed energy En = 789 cm^^ is 3% lower than the 
experimental value =811 cm^^, whereas for n = 35 
our computed energy En = 889 cm^^ is 5% lower than 
the experimental value En ~ 940 cm~^. This divergence 
in the upper part of the spectrum is undoubtedly cor- 
related to the fact that our dissociation energy is lower 
than the experimental dissociation energy. 

To see how the normalized energies change with respect 
to variations in [i and a, we recompute the potential en- 
ergy curve, as well as the vibrational spectrum supported 
by it, for various values of [i and a in the neighborhood 
of /i = 1/1836 and a = 1/137. The computational chem- 
istry programs assume atomic units; numerical variations 
in /i and a are effected by modifying parameter values 
for the reduced mass within our Matlab routine and the 
speed of light within DIRACIO [M = Ufj.'^ a.u. for 
^*Sr2 and c= 1/a a.u., where a.u. denotes the respective 
atomic units of mass and velocity) . We then obtain the 
sensitivities dfj,£n and da£n from numerical differentia- 
tion with respect to fi and a. We emphasize that our 
method for obtaining these sensitivities treats variations 
in fi and a in a similar manner and on equal footing. Fig- 
ure [3] displays our results for d^^Sn and daSn for the levels 
n = through n = 35; these level sensitivities arc plot- 
ted versus the normalized energy £„ . We note a behavior 
which resembles that "predicted" by the Morse potential. 
Namely, both sensitivity curves approach the appropri- 
ate limits for — >■ and f„ — >■ 1, while simultaneously 
peaking at = 3/4. Moreover, the curves are found 
to be essentially proportional, with the ratio daSn/d^Sn 
being 0.28±0.02 across the entire range of data (and 0.28 
at the common maximum of the two curves). 



V. CONCLUSION 

Here we have considered the influence of a variation 
on the experiment proposed by Zelevinsky et al. 0, Q to 
probe variation in the electron-to-proton mass ratio using 
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FIG. 3. (Color online) Level sensitivities d^iEn and da£-n vs 
fn for the X state of Sr2. The circles correspond to ah 
initio results for individual levels n = through n = 35, with 
the solid curves being fits to this data (the fits are extended 
through £n = 1). 



the vibrational spectrum of Sr2 . The relevant observable 
in this experiment is the frequency ratio R, illustrated in 
Fig. [TJ With one "anchor" level taken at the bottom of 
the spectrum and another at the top, the frequency ratio 
is given approximately by 

R ^ 2£n ~ 1 , 

where n labels the intermediate level (e.g., n — 27 in 
Fig. [1}. Our ab initio computations predict that the 
frequency ratio R is only 3 to 4 times less sensitive to 
variation in a as it is to variation in /i. Specifically, we 
find that variations in /i and a induce a variation in R 
according to the relation 

SR^Kf^Ji +0.28^-^). (3) 

We estimate the uncertainty in the factor 0.28 to be on 
the order of 10%, based primarily on the discrepancy 
of our computed dissociation energy to the experimental 
dissociation energy. Equation ^ summarizes the prin- 
ciple result of this work. The factor K here is given ap- 
proximately hy K = 2 X df^Sn] non-zero shifts in the two 



anchor levels amount to small corrections which reduce 
K from this value. 

As suggested by Eq. (jS]), a measured drift in the fre- 
quency ratio R cannot, by itself, be used to distinguish 
between fi variation or a variation. To extract informa- 
tion about variations in the respective constants them- 
selves requires further experimental input. Optical ion 
clocks have been used to test the stability of a, with the 
ratio of clock frequencies being insensitive to ^ variation. 
The current best limit on a variation allows for a drift 
on the fractional level of 4 x 10~^''/year For the pro- 
posed experiment of Zelevinsky et al., this result may be 
used with Eq. ([3]) to justify neglect of a variationjwhich 
has been implicitly assumed in previous works 0, Q • On 
the other hand, for high experimental precision — namely, 
experimental precision capable of detecting a drift in fi 
at the fractional level of 1 x IG^^^/year — equation ([3]) 
indicates that a variation should not be neglected. Such 
high precision is conceivable; in the related proposal of 
DcMillc et al. [1], referred to in the Introduction, the 
authors argued that their method could plausibly detect 
fractional variations in /x at < 10^^^. 

For such high experimental precision, additional exper- 
imental input could perhaps be obtained by substituting 
Sr2 with another species, such as Yb2, in the experiment. 
Yb has a similar valence structure as Sr and also has 
isotopes which lack nuclear spin (168,i70,i72,174,i76y13)^ 
Moreover, as with Sr, high precision spectroscopy on 
optically trapped Yb has become a refined art [l^, [l3 ■ 
From the {aZY scaling of the relativistic corrections to 
the electronic potential, we may presume that an Yb2 
experiment would be about equally sensitive to a varia- 
tion as to /i variation, with an estimated sensitivity ratio 
0.28 X (70/38)2 ^ Q \]s\n.g Sr2 and Yb2 results in con- 
junction, one could conceivably determine both /x and a 
variation to high precision with the proposed experiment 
of Zelevinsky et al. 
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